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Abstract: We study the problem of computing the rank of a divisor on a finite graph, a quantity 
that arises in the Riemann-Roch theory on a finite graph developed by Baker and Norine {Advances 
in Mathematics, 215(2): 766-788, 2007). Our work consists of two parts: the first part is an algorithm 
whose running time is polynomial for a multigraph with a fixed number of vertices. More precisely, our 
I—, algorithm has running time 0(2"'°s")poly(size)(G), where n + 1 is the number of vertices of the graph 

G. The second part consists of a new proof of the fact that testing if rank of a divisor is non-negative 
or not is in the complexity class NP Hco— NP and motivated by this proof and its generalisations, we 
• construct a new graph invariant that we call the critical automorphism group of the graph. 

a 

^ 1 Introduction 

The study of chip firing games on a graph with their several variants is a classic topic in discrete 
^ mathematics, in particular in algebraic graph theory and this is reflected in the fact that standard text 

books on algebraic graph theory, such as the book by Godsil and Royle j7] have chapters devoted to 
chip firing games. See the article of Merino ^14] for a broad survey of the topic. 

Despite their simplicity, chip firing games have connections to various other areas of mathematics 
' and physics. Some examples include chip firing games studied in dynamical systems under the name 

"sandpile models" and in fact, early algorithmic work on chip firing games was done independently by 
y—i physicists. Chip firing games of graphs also play a role in counting the number of points on elliptic 

curves over finite fields and is the topic of the Phd thesis of Greg Musikar. Recently, chip firing games 
^ have played a key role in developing analogies between graphs and Riemann surfaces. This line of 

k>J research was initiated by Baker and Norine in their pioneering work "Riemann-Roch and Abel-Jacobi 

rS theory on a finite graph" |2] where they show analogues of the Riemann-Roch theorem on a finite graph 

^ and the theorem is best explained in the language of chip-firing games. 

A chip-firing game is a solitary game played on an undirected connected multigraph and is defined 
as follows: Each vertex of the graph is assigned an integer, refereed to as "chips" and this assignment 
is called the initial configuration. At each move of the game, an arbitrary vertex v is allowed to either 
lend or borrow one chip along each edge incident with it and we obtain a new configuration. We define 
two configurations Ci and C2 to be equivalent if Ci can be reached from C2 by a sequence of chip 
firings. The Laplacian matrix Q of the graph naturally comes into the picture as follows: 

Lemma 1. Configurations Ci and C2 are equivalent if and only if Ci — C2 can he expressed as Q ■ w 
for some vector w with integer coordinates. 



We can ask some natural questions on such a game: 



1. Is a given configuration equivalent to an effective configuration i.e., a configuration where each 
vertex has a non-negative number of chips? 

2. More generally, given a configuration what is the minimum number of chips that must be removed 
from the system so that the resulting configuration is not equivalent to an effective configuration? 
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The Riemann-Roch theorem of Baker and Norine provides insights into answering these questions. 
We need the following definitions before we can state the theorem. 

Definition 1. (TDivisor^ A configuration on a finite graph is called a divisor and is represented as an 
integer vector on n + 1 coordinates. 

Definition 2. (Degree of a divisorj For a divisor D, the total number of chips, i.e., the sum of chips 
over all the vertices of the graph, is called the degree of the divisor D, denoted by deg{D). 

Definition 3. (Hank of a divisor j For a divisor D, one less than the minimum number of chips that 
must be removed from the divisor D so that the resulting configuration is not equivalent to an effective 
divisor is called the rank of the divisor D, denoted by r{D). 

Theorem 1. (lliemann-Roch theorem for graphs^ For any undirected connected graph G, there 
exists a divisor K called the canonical divisor such that for any divisor D we have the following formula: 

r{D) - r{K - D) ^ deg{D) - {g - 1) (1) 

where g is the cyclotomic number of G and is equal to ni — n + \ where m is the number of edges, n is 
the number of vertices of G. 



1.1 Related work and a brief description of our work 

The central quantity in the Riemann-Roch theorem is the rank of a divisor and the efficient computation 
of rank is a natural problem, attributed to Hendrik Lenstra (see [3]). In fact this work, gave an algorithm, 
i.e., a procedure that terminates in a finite number of steps to compute the rank of a divisor on a tropical 
curve. But the algorithm docs not run in polynomial time in the size of the multigraph even when the 
number of vertices are fixed since the algorithm involves iterating over all the spanning trees in the 
graph (see Proof of Theorem 23 in [4 ) and the number of spanning trees in indeed not polynomially 
bounded in the size of the mutligraph even if the number of vertices are fixed. On the other hand, 
there are polynomial time algorithms for deciding if the rank of a divisor on a finite multigraph is 
non-negative, see [6], [17] and [15]. 

Our paper is centered around the problem of computing the rank. More precisely, (Section [s]) we 
obtain an algorithm whose running time is polynomial for a multigraph with a fixed number of vertices. 
More precisely, our algorithm has running time 2'^^"'°sn)poly(size)(G), where n -f 1 is the number 
of vertices of the multigraph G. Recall that we are working with arbitrary undirected connected 
multigraphs or equivalently graphs with positive integer weights on the edges and indeed, the original 
Riemann-Roch theory was also developed in this setting. The main tools involved are the Riemann- 
Roch formula and a formula for rank (Theorem [4]) that is used in the proof of the Riemann-Roch 
formula, these results were first obtained in the work of Baker and Norine [5]. We obtain a geometric 
interpretation of rank (Theorem [s]) and combine this geometric interpretation along with algorithms 



from the geometry of numbers to obtain the algorithm for computing the rank (Algorithm 3.5 ). We find 
it satisfying that geometric tools seem to be essential in obtaining the algorithm, though the definition 
of rank of a divisor can be stated in purely combinatorial terms. 

The second part (Section |4]) starts with a new proof of the fact that testing if r{D) > is in 
NP n CO — NP. A duality theorem characterising divisors with r{D) > plays a key role in our 
proof. Motivated by the observation that generalisations of the duality theorem lead to more general 
complexity results on computing the rank, we generalise the duality theorem and this generalisation 
leads to the construction of a new graph invariant that we call the critical automorphism group of 
the graph. 



2 Preliminaries 



We will now describe mainly geometric notions and results that we frequently use in the rest of our 
paper. 
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2.1 Lattices 



A lattice i is a discrete subgroup of the Euclidean vector space E". More concretely, a lattice is the 
Abelian group obtained by taking all the integral combinations of a set of linearly independent vectors 
61 , . . . , 6fe in ]R" . More precisely, 

fc 

L^{J2aA\ (2) 

i=l 

The set B = {61, . . . , bk} is called a basis of L and the integer k is called the dimension of L, denoted 
by dim{L), is independent of the choice of the basis. We denote the subspace spanned by the elements 
of B by Span{L). 

Let us now look at the important geometric invariants of a lattice. The volume of the lattice Vol{L)^ 
also known as the discriminant or determinant, is defined as ^/ det{BB*) where the basis B is represented 
as a matrix with its elements row-wise and hence, BB* is the Gram matrix of the basis elements. Another 
important invariant is the norm of a shortest vector of a lattice. 

Definition 4. A shortest vector of L in the Euclidean norm is an element q of L such that q-q< q' -q' 
for all non-zero elements q' in L and we denote \\q\\2 as ve{L). 



2.2 The Laplacian Lattice of a Graph 



For an undirected connected graph G, the Laplacian matrix Q{G) is defined as D(G) — A{G) where D{G) 
is the diagonal matrix with the degree of every vertex in its diagonal and A(G) is the vertex-adjacency 
matrix of the graph. We assume the following standard form of the Laplacian matrix: 



So 

-bio 



-boi 
Si 



-bo2 ■ 
-bi2 . 



— bnO ^bnl —bn2 



-bon 
-bin 

Sn 



(3) 



has the following properties: 

(Ci) bij^s are integers, bij > for all < i =/= j < n and bij ~ bji, Vi =/= j. 
{G2) Si = YTj=i,j=^i bij = YJj^i.j^i b]i (and is the degree of the i-th vertex). 

Though the Laplacian matrix of a graph contains essentially the same information as the adjacency 
matrix of a graph, it enjoys other nice properties. See, for example, the chapter "The Laplacian matrix 
of a Graph" in the algebraic graph theory book of Godsil and Royle i7j for a more complete discussion. 

Lemma 2. The Laplacian matrix Q{G) is a symmetric positive semi-definite matrix. 

Another remarkable property of the Laplacian matrix is described in the Matrix- Tree theorem: 

Theorem 2. (KirchofF's Matrix- Tree Theorem) The absolute of value of any cofactor is equal to 
the number of spanning trees of the graph. 



A remarkable aspect of the matrix-tree theorem is that it reduces counting the number of spanning 
trees of a graph into a determinant computation and hence provides a polynomial time algorithm for 
it. 

Definition 5. (The Laplacian lattice of a graph) Given the Laplacian matrix Q{G), the lattice 
generated by the rows (or equivalently the columns) of Q{G) is called the Laplacian lattice Lq of the 
graph. 

In the language of chip firing games, the Laplacian lattice is the set of all divisors that are equivalent 
to the divisor (0, . . . , 0). 



3 



Definition 6. (The Hyperplanes Hk) For a fixed real number k, we denote n- dimensional hyperplane 

{{xo,...,xn)\ Er=o^« = ^} by Hk- 

Definition 7. (The Root Lattice yl„) The root lattice is the lattice of integer points in the 
hyperplane Hq = {{xq, . . . , a;„)| J^i^o — G M} = (1, . . . , 1)^. More precisely, 



Remark 1. The name "root lattice'' is derived from the fact that the lattice An is generated by a root 
system i.e. a set of vectors that satisfy reflection symmetries. See pages 96-98 of Conway and Sloane 
for a precise definition of a root system and for a discussion on root lattices. In the case of An, the 
corresponding root system is Ci — ej where Ci and Cj run over the standard basis ofW^^^. 

We now make a few simple observations on the Laplacian lattice of a graph. 

Lemma 3. The Laplacian lattice of a graph on n + 1-vertices is a sublattice of the root lattice An. 

Definition 8. (The covolume of a sublattice) A full- dimensional sublattice Lg of a lattice L is a 
subgroup of the Abelian group L the cardinality of the quotient group L/Ls is called the covolume of Ls 
with respect to L. 

Lemma 4. /J 5/ The covolume of the Laplacian lattice of G with respect to An is equal to the number 
of spanning trees of G. 

The elements of An/ Lq naturally possess an Abelian group structure and this group is known as 
the Picard group of G, also known as the Jacobian of G. A number of works have been devoted to the 
study of the structure of this group and the information that it contains about the underlying graph, 
see for example the works of Biggs [5] , Kotani and Sunada |TT] and Lorenzini [T? . As a straightforward 
corollary to Lemma |4] we obtain: 

Corollary 1. The cardinality of the Picard group of G is equal to the number of spanning trees ofG. 
2.3 Algorithmic Geometry of Numbers 

We will now briefly discuss some important algorithmic problems related to lattices and geometry of 
numbers, in particular those that we employ in our algorithm. Two important and closely related 
problems in the algorithmic geometry of numbers are: 

• Closest vector problem (CVP): Given a lattice L by an arbitrary basis and a target vector 
V, find a lattice that is closer to v than to any other lattice point in the i'2-norm. 

• Shortest vector problem (SVP): Given a lattice L, by an arbitrary basis find a non-zero 
lattice point that has the smallest ^2-iiorm. 

Indeed, CVP and SVP have versions with respect to the other £p-norms. For the ^2-iiorm both CVP 
and SVP are known to be NP-hard and in fact SVP has a polynomial time reduction to CVP in any 
given norm but the converse is not known. A problem that is closely related to CVP is the integer 
programming problem: 

Integer programming problem: Given a polyhedron V in the form A-x <b where A is an m x n 
matrix and 6 is a vector in M'". Decide if V has an integer point or not i.e., 7-" n Z" = or not. 

Indeed, the integer programming problem is also known to be NP-hard. In the 1980's there was a 
great amount of progress in algorithmic geometry of numbers, trigerred by the algorithm of Lenstra 
that solves the integer programming problem in polynomial time for a fixed dimension [TU]. Lenstra's 
algorithm has running time 2*^'" • poly{\I\), where |/| is the size of the input. The factor 2'^'" ^ 
was subsequently improved and the current best factor being 2'^("'°sn) by Kannan [S]. Note that the 
Kannan's algorithm also works when the polytope is presented as a separation oracle and in fact, the 
algorithm works also for general convex bodies presented in terms of a separation oracle, see the remark 
"General convex bodies and mixed integer programs" in f8|. We will crucially use Kannan's algorithm 
and the polytope will be presented to Kannan's algorithm in terms of a separation oracle. 



n+l 




4 



2.4 Polyhedral Distance Functions 



Let V he a. convex polytope in E" with the reference point O = (0, . . . , 0) that we cah the "center" in 
its interior. By V{p, A) we denote a dilation of P by a factor A and its center translated to the point 
p i.e. V{p, X) = p + X.V and X.V = { X.x \ x & V}. We define the P-midpoint of two points p and q 
in M" as inf{i?| 'P{p, R) fl V{q, R) ^ 0}. The polyhedral distance function d-p{. , .) between the points of 
M" is defined as follows: 



d-p is not generally symmetric, indeed it is easy to check that d-p(. , .) is symmetric if and only if the 
polyhedron V is centrally symmetric i.e. V = —V. Nevertheless d-p{. , .) satisfies the triangle inequality. 

Lemma 5. For every three points p,q,r G M", we have d-p{p,q) + d-p{q,r) > d-p{p,r). In addition, if 
q is a convex combination of p and r, then d-p{p,q) + d-p{q,r) = d-p(p,r). 

Proof. To prove the triangle inequality, it will be sufficient to show that ii q E p + X.V and r E q + 1.1.V, 
then r E p + {X + i.j,).V. We write q = p + X.q' and r = q + fi.r' for two points q' and r' in V. We can 
then write r = p + X.q' + ^.r' = p + {X + m)( aTa^ ^ a+TI'^')' ^ being convex and A, ^ > 0, we infer 
that yr— .g' + Y^-r' E V, and so r E p + (X + 11). V. The triangle inequality follows. 

To prove the second part of the lemma, let t E [0, 1] be such that q = t.p + (1 — t).r . By the triangle 
inequality, it will be enough to prove that d-p{p,q) +d-p{q,r) < d-p{p,r). Let d-p{p,r) = A so that r = 
p+X.r' for some point r' in V. We infer first that q = t.p+{l—t).r = t.p+{\—t){p~\-X.r') = p+{l—t)X.r\ 
which implies that d-p{p, (?) < (1 — t)X. Similarly we have t.r = t.p + tX.r' = q — {1 — t)r + tX.r' . It 
follows that r = q + tXr' and so d-p{q,r) < tX . We conclude that d-p {p, q) + d-p{q,r) < d-p {p,r), and the 
lemma follows. □ 

We also observe that the polyhedral metric d-p{. , .) is translation invariant, i.e. 

Lemma 6. For any two points p, q in M", and for any vector v E M", we have d-p(p, q) = d-p{p—v, q—v). 
In particular, d-p {p, q) = d-p{p — q,0) = d-p{0,q — p). 

Proof. The proof is easy: \i q E p + X.V , then q — vEp — v-{- X.V, and vice versa. □ 

Remark 2. The notion of a polyhedral distance function is essentially the concept of a gauge function 
of a convex body that has been studied in [16] . Lemmas [s] and [6] can be derived in a straight forward 
way from the results in |16| . 

Recall that the n-dimensional hyperplane i?fc is defined as Hi^ = {(xg, . . . ,Xn)\ Yl^=o^i — 
will be mainly be using distance functions defined by regular simplices A and A where the simplices A 
and A are defined as follows: 

Deflnition 9. The regular simplex A is the convex hull of to, . ■ . ,tn E Hq where 



Vp,geM", drip,q) : 



M{X>0\qEVip, A)}. 




1, otherwise 



for i from 0, . . . ,n and tij is the j-th coordinate ofti. We define A as —A. 



We note that for points in Hq, the distance functions d^ and have a simple formula: 



Lemma 7. (Lemma 4-7, 111) For any pair of points p, q in Hq, we have: 



dA{p,q) = |min((7j - p^)\, 



dA{p,q) = |min(K - qi)\. 
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Given a permutation tt on the n + l vertices of G, define the ordering 7r(i;o) <Tr 7r(wi) • • • <7r T^ivn) 
and orient the edges of graph G according to the ordering defined by tt i.e., there is an oriented edge 
from Vi to Vj if (vi, Vj) € E and if Vi <,r Vj in the ordering defined by tt. Consider the acyclic orientation 
induced by a permutation tt on the set of vertices of G and define = {indegT^{vQ) — l, . . . , indegT^ivn) — 
1), where indegjriv) is the indegree of the vertex v in the directed graph oriented according to tt. Define 

Ext(LG) = {-i^TT + q\TT e Sn+i, q G La}- (4) 

2.5 Distance function induced by a Discrete Point Set 

Given a polyhedral distance function V and a discrete point set S, we define a function h-p s ■ — > R 
as: 

hv,s{p) = mm{d-p{p, q)} (5) 

q£S 

In particular, the notion of local minima and local maxima of the distance function h-p,s turns out 
to be useful: 

Definition 10. (Local Maxima and Local Minima of h-p s) Let B{p,e) be the Euclidean ball of 
radius e centered at p. A point c in M" is called a local minimum of h-p_s if there exists an e > such 
that h-p^s{c) < h-p^siq) for all q € B(c,e). A point c in M" is called a local maximum of h-p, s if there 
exists an e > such that h-p^s{c) > h-p^siq) for all q e i?(c, e). 

We denote the set of local maxima of h-p^s by Crit7:>(5). 

3 Algorithms for computing the rank 

In this section, we will construct algorithms for computing the rank with the main result being an 
algorithm for computing the rank that runs in polynomial time when the number of vertices of the 
multigraph is fixed. 

3.1 A simplification 

We shall first observe that by using the Riemann-Roch theorem we can restrict our attention to divisors 
of degree between zero and g ~ 1- Firstly, a divisor of negative degree must have rank minus one. 
Furthermore, by the Riemann-Roch formula we have: 

Lemma 8. // the degree of D is strictly greater than 2g — 2, then r(D) — deg{D) — g. 

Proof. Observe that if the degree of D is strictly greater than 2g — 2 then the rank of K — D is —1 and 
apply the Riemann-Roch theorem. □ 

Furthermore, we can compute the rank of divisors of degree between g and 25 — 2 by computing the 
rank of K — a divisor that has degree between zero and .g — 1 and then applying the Riemann-Roch 
theorem. Hence, we consider the problem of computing the rank of a divisor of degree between zero 
and 5 — 1. In fact, we consider the decision version of the problem i.e., we want to decide (efficiently) 
if r{D) < k for every k between zero and 5 — 1; observe that such a procedure combined with a binary 
search over the parameter k will compute the rank in time 0{h\{g)) times the running time of the 
procedure. 

3.2 A first attempt at computing the rank 

Let us discuss a first attempt at computing the rank. We will compute rank directly from its definition 
(Definition [3]) . We will use the fact that there is a polynomial time algorithm for testing if r{D) > 
due to the independent work of Dhar and Tardos [TT] . 
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Algorithm 1. 1. Enumerate all effective divisors of degree at most the degree of the divisor D . 

2. Find an effective divisor E of smallest degree such that r{D — E) = —1 by using Dhar's algorithm. 
Theorem 3. The running time of Algonthm\^is 0(2" i"^). 

The running time of the Algorithm [T] is not polynomial in the size of the input even for a fixed 
number of vertices since the quantity 2"'"^* is not polynomially bounded in the size of the input. 
There is general interest in obtaining an algorithm that runs in polynomial time for a fixed number 
of vertices and furthermore, in obtaining a singly exponential time algorithm i.e., an algorithm with 
running time 2'^'^"'poly(size(G)). We will now undertake a deeper study of rank to obtain an algorithm 
that runs in time polynomial in the size of the input provided that the number of vertices is fixed. 
More precisely, our algorithm has running time 2'^'^"'°s")po/?/(size(G)). An important ingredient is a 
geometric interpretation of rank that we shall obtain in the following section. 

3.3 A geometric interpretation of rank 

We start with the following formula for rank first shown in Baker and Norine [2^ and later reproven in 
Amini and Manjunath [T]. 

Theorem 4. For any divisor D, we have: 

r{D) ^ min^fzExt{LG)deg^{D - ly) - 1 (6) 

where deg+{D) = J2i:D,>o^i- 

Remark 3. Some remarks on the proof (s) of Theorem^are in order: as we mentioned earlier, Theorem 
has two proofs, the original proof due to Baker and Norine was based on combinatorial tools. In 
particular, the main component of the proof was to establish the existence and uniqueness of a certain 
special type of divisors called "v-reduced" divisors in each linear equivalence class of divisors, while 
the approach of Amini and Manjunath IjJ involved studying the Laplacian lattice under the simplicial 
distance function d/x ■ 



3.3.1 A sketch of the approach 

Let us now briefly sketch our approach to computing the rank: We start with the formula to compute 
the rank and proceed as follows: we run over all the permutations tt G Sn+i and for each permutation 
TT suppose that we could compute miug^j^^ deg^{D — f tt + in time that is possibly exponential but 
only in n then we would obtain an algorithm with running time 0{f{n)poly{size{G))) for some function 
/. But, how do we compute miuggi^ deg^{D — + q)l One hope would be to reduce the problem to 
a closest vector problem on lattices or more generally to integer programming. Fortunately, the integer 
programming problem has an algorithm that runs in time that exponential only in n (the dimension of 
the lattice). Such an algorithm would run in time 0(2"i°8>o;?/(size(G))). This approach requires a 
better understanding of the deg^ function that we now obtain. 

Definition 11. (Orthogonal projections onto Hk) For a point P G ]R"+^ we denote by iTk{P) the 
orthogonal projection of P onto the hyperplane Hk . 

For the sake of presentation, we first consider the case where the divisor has degree g — 1. In this 
case, we observe that deg^{D — v) = ^ and that £i{D — v) = £i(7ro(D) — 7ro(i^)). We denote 

the set TTo{Ext{Lc)) the orthogonal projection of Ext{Lc) onto the hyperplane Hq by CritA(iG), as 
defined in Subsection |2.5| and indeed, the orthogonal projections of Ext{LG) are the local maxima of 
the distance function /ia.Lg refer to [1] for more details. 

Corollary 2. For any divisor D with deg{D) = g — 1, we have: 

r{D) = mm 1. (7) 

ceCritA(LG) ^ 
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Figure 1: The distance function 



Taking cue from Corollary [2] it is natural to ask if there is a similar "distance function" type 
interpretation for divisors of degree between zero and g — 1. We will answer this question in the 
affirmative, the relevant distance function, actually a family of distance functions is the following: 

Definition 12. (Degree-Plus Distance^ Let fc be a positive real number. For points P and Q in 
Hq, we define the generalised degree-plus distance between P and Q as 

=sup{r| A(P,r)n A(Q,r + fc) =0}. 
Note that though df does not appear to be a distance function at first glance, we will actually show 



that it is can be realised by a sequence of distance functions (See Section 2.4 for a definition) 
We will now note some basic properties of the function d^: 

Lemma 9. (Translation InvarianceJ For any points P, Q, and T in Hq and for any positive real 
numbers ri and r2 we have: A{P, vi) C] A{Q, r2) — % if and only if A{P + T, ri) n A{Q -\- T, r2) — 0. 

Proof. Assume that A(P,r-i) n A{Q,r2) ^ and consider a point S E A(F, ri) n A(Q,r2). Now, S = 
ri Yli=i Oii'^i + P = r2 Yh=i Pi'^i + Q for some > 0, /3i > and = Si ft = 1- Now this imphes 

that S+T = n ^"+/ a,v,+P+T = ra Er=i' P^v,+Q+T. Hence, e A{P+T, ri)nA{Q+T, ra) ^ 0. 
The converse follows by symmetry. □ 

Lemma 10. (Projection Lemmaj Let R be a point in IR"+^ with deg{R) > 0, let O be the origin 
and let 7To{R) be the orthogonal projection of R onto Hq. We have: 

inf deg(Z) = (n + 1) sup ( r I A(7ro(i?), r) n A(0, r + ^^f^i^) = 0\ + deg(R). 

ZeH+{R)nH+{0) J f ^ t y ^, , v , „ _^ 1 ^ J 

Proof. Consider a point X, say, in the intersection of H^{0) and H^{R) and consider the intersection 
of the hyperplane H(if,g{x) with H^{0) and H^{R). Observe that the intersection of H^pg^^x) and 
H^{0) is a simplex that is a scaled and translated copy of A, call it Ai, centered at ^^^^^^(1, ■ ■ ■ , 1) 
and scaled by a factor of ^^^^^^3^. Similarly, the intersection of -ff(ieg(x) and H^{R) is also a simplex that 
is a scaled and translated copy of A, call it A 2 centered aX R+ ■ • • , 1) scaled by a factor of 

''^^1+]"^^ . Observe that deg{X) > deg{R) > deg{0). Indeed simplices Ai and A2 intersect at X and 
deg{X) is equal to n + 1 times the radius of A2 plus deg{R). We now project the simplices Ai and A2 
onto ITo and obtain inf {de5(Z) I Z e H+{R)nH+{0)} > (n + l)sup{r| A(7ro(i?), r)n A(0, r+ ^fi^) = 
0} + deg{R). Now, consider a point P in A(7ro(i?),r) n A{0,r + ^^^^pj^) and observe that the point 
X = P + {r + . . . , 1) is a point in the intersection of H+{0) and H+{R). This shows that 

mi{deg{Z)\ Z e H+{R)nH+{0)} < (n + l)sup{r| A{7To{R),r)n A{0,r + ^^^) ^ (!)} + deg{R). This 
completes the proof. □ 



We are now ready to establish the connection between the deg^ function and the function d 



k ■ 
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Lemma 11. For any pair of points P and Q in M."''^^ with deg{P) > deg{Q), we have 
deg+iP^Q) = in + l)d+{7ToiP),MQ)) + deg{P-Q) 

for k = ^-Ml^ 



Proof. First consider a point R in R"-+^ We have deg+{R) = J2r,>o = deg{R®0), where R®0 
{max{Ro, 0), . . . , max{Rn, 0)). Now we have: 

deg{R®0)= inf deg(Z). 

zeH+(R)nH+(o) 

By Lemma [To] we have 

inf degiZ) = (n + 1) sup ( r I A(7ro(i?), r) n A(0, r + ^^^^) = 0| + deg(i?). 

ZeH+{R)nH+(0) J f ^ I \ vy ; \ , „ + l ' J ^ 



Now, for two points P and Q in ]R"+^, letting i? = P — Q in the above formula and applying Lemma 
[9| we obtain the relation given in the proposition. □ 

The function is motivated naturally by the definition of the deg'^ function but is not very handy 
for geometric as well as computational reasons. In the following, we will obtain a more convenient 
representation of d^. In fact, is closely related to the following family of polytopes: For a point 
P G Hq and m,n > 0, let 'Pm,n{P) = (A(0,m) © 

Mink ^{0,Ti)) + -P, whcrc ®Mink denotcs the 
Minkowski sum. Note that we use the notation © for the tropical maximum sum. 

Lemma 12. For any positive real numbers m, n, Pm,niP) « convex polytope. 

Proof. Using the fact that Minkowski sum of two convex polytopes is a convex polytope and hence, 
'Pm,n{0) is a convex polytope. Indeed translates of a convex polytope is also a convex polytope and 
hence, Vm,n{P) is also a convex polytope. □ 

Lemma 13. For points P and Q in Hq and for k>0, d'^{P, Q) = inf{r| Q e (A(0, r) ® Mink ^{O ,r + 
k)) + P}, where O is the origin. 

Proof By definition d^ (P, Q) = sup{r\ A(P, r) n A(Q, r + k) = 0}. Let ro = d^ (P, Q) and consider a 
point R in the intersection of A(P, rp) and A{Q,ro + k). Rephrasing dA(P, i?) — ro and d/\{Q,R) = 
(P, Q) = ro + k. This imphes that P - P G A(0, ro) and Q - R e A(0, ro + fc). This shows that 
Q - P e A(0, ro) ©Aftnfc A(0, ro + fc) and we obtain Q e (A(0, ro) ©Mmfc A(0, ro + fc)) + P. Hence, 
d+l^^, Q) > inf{r| g e (A(0, r) ©Mmfc A(0, r + fc)) + P}. 

Furthermore, if Q is contained in (A(0,r) ©A/mfe A(0,r + fc)) + P then there exists a point P = 
Pi + P2 such that Pi e A(0, r) and P2 e A(0, r + fc) with Q = Pi + P2 + P and we take P3 such 
that P3 = Q — P2 = Pi + P. Therefore, the point P3 is contained in both A(Q, r + fc) and A(P, r) and 
we obtain inf{r| Q e (A(0, r) ©Mmfe A(0, r + fc)) + P} > d;!:(P, Q). This concludes the proof. 

□ 



As a corollary we obtain a handy characterisation of the polytope 'Pm,n' 

Corollary 3. Let P, Q be points in Hq, a point Q belongs to the polytope Vr,r+d{P) if o,nd only if 
deg+{P+'^^^Q)<r{n + l). 

Proof. Let a point Q be contained in 'Pr,r+d{P), then inf{r'| Q e {/S{0,r')®Mink^{0, r' + d)) + P} < r 
and hence by Lemma [Ts] we know that r > inf{r'| Q e (A(0, r') ©M^nfe A(0, r' + d)) + P} = d+{P,Q). 



ByLemmajn] weknowthat deg+(P+^^%^-0) = (n+l)dJ(P, Q). Hence, deg+(P+%^-Q) < 
r{n + 1). Conversely, if a point Q satisfies deg^{P + '^''^'^i"'"'* ~ Q) < r{n + 1) then, d^{P, Q) < r and 
hence, r > inf{r'| Q £ (A(0,r') ©Mmfc A(0, r' + d)) + P} and hence Q £ Vr.r+d{P). □ 
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Now putting together, the formula for rank in Theorem [4j Lemma 11 and Lemma 13 we obtain the 
foUowing geometric interpretation of rank: 

Theorem 5. (A Geometric Interpretation of rankj Consider a divisor D of degree d between zero 
and g ~ I, then D has rank tq — 1 if and only if ttq{D) is contained in the boundary of the arrangement 
Uc6CritA(LG)'^ri,r2(c) where ri = ro/{n + 1) and r2 = (ro + .g - 1 - d)/{n + 1). 

Remark 4. The fact that n+kic) C r^+kic) if ri < r2 is implicit in the statement of Theorem 

m 

3.4 Computing the rank for divisors of degree between zero and g — 1 

We now give an algorithm for computing the rank that runs in polynomial time for a fixed number of 
vertices. The algorithm uses two main ingredients: 

1. The geometric interpretation of rank (Theorem [5|. 

2. Reduction to the algorithm for integer programming by Kannan [5]: 

For the sake of exposition, we first consider the slightly easier case of divisors with degree exactly 
g — 1. We employ Theorem [2] to obtain the following algorithm: 

Algorithm 2. 1. For each permutation tt G Sn+i, we compute min^gi^ £i(ttq(D) — q)/2 — 1 using 
Kannan's algorithm. The ii unit ball is given to Kannan's algorithm as a separation oracle (we 



will provide an efficient implementation of the separation oracle in Lemma 14 )■ 
2. We minimise over all permutations tt. 

We now turn to the general case: we start with the geometric interpretation for rank and we would 
like to reduce the problem to the integer programming problem. We construct a preliminary algorithm 
as follows: 

Algorithm 3. 1. Find the smallest integer r such that 7ro(I?) is contained in UcgCritA(iG)^»'i>''2 (c) 
where ri — , r2 = ' ^^^+1^^^ ^csim;; for all values of r from zero to g ^ 1. 

Using the fact that the degree of the divisor is between zero and g — 1, the algorithm would run in 
time 0{g ■ 2'^("^°s") • poly{size{G))). Since g is not polynomially bounded in the size of the input, the 
algorithm does not run in polynomial time for fixed values of n. We resolve this problem by performing 
a binary search over the parameter r in the polytope Vri.r^i^^) ^-nd s-Pply Kannan's algorithm at each 
step of the binary search. Since we know from Subsection |3.1[ that the rank of the divisor is at most 

1 the algorithm terminates in 0(2"'°s"poZy(size(G'))). Here is a formal description of the algorithm: 

Algorithm 4. 1. For each permutation n e Sn+i, use binary search on the parameter r along with 
Kannan's algorithm to test if 'Kq{D) is contained in Uq^L^Vr^r+g-i-dicrr + q), the polytope is 
presented to Kannan's algorithm as a separation oracle. 

2. Repeat over all permutations tt. 

The straightforward way of presenting the polytope Vr.r+k to Kannan's algorithm is in terms of 
its facets. But since the number of facets of Pr,r+k is 2"+^, the factor depending on n in the time 
complexity of the algorithm becomes larger than 2"'°s". Hence, we present the polytope Vr.r+k by a 
separation oracle to Kannan's algorithm and the following efficient implementation of the separation 
oracle ensures that the algorithm runs in time 2'-"-"^°^^^^poly{size{G)). 

Lemma 14. (A separation oracle for the polytope Vm^n) There is a polynomial time separation 
oracle for the polytope Vr.r+d given any point p there is a polynomial time (in the bit length of 
p and the vertex description of Vr,r+d) algorithm that either decides that p is contained in Vr,r+d or 
outputs a hyperplane separating the point p and the polytope Vr,r+d- 
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Proof. Given a point Q, compute the function r' = deg~^{P + '^^n^l'i'^'' ~ Q) /("^ + 1) a-nd if r' < r then Q 
is contained in Pr,r+d{P) and otherwise let 5+ be the set of indices such that Pi + d{l, . . . , l)/n + 1 > 
0, output the hyperplane Hs'- J2x GS+(~^i + Pi + d{l, . . . ,1) /n + 1) < ('' +'"K"+i) separating 
hyperplane. To show that Hg is a separating hyperplane, assume the contrary and since X^a- gs+ ("•^j ~^ 
P^+^^^^) = r'{n+l) > there isapointQin7^^^„ such that XQ^gg+(P, + ^'^j^-gO > 

(«+iKr'+r) ^ ^j^^^ p^j^^ qL ^ Q _ ^ih^ we have: Ejes+(<3^)j = ^'(" + 1) > 

^ jj^^^^^ deg+{Q^) > ('•+'-'K"+i) > r(n + 1) and deg(g^) = d, using Corollary |3| we obtain 
a contradiction. 

□ 

The correctness of the algorithm is clear from Theorem [sj 

Theorem 6. For any divisor D with degree between zero to g ~ 1, Algorithm^ computes the rank of 
the divisor D. 

Theorem 7. ^/(7oni/i77i[^ runs in time 2'~'^''^^°^'^^poly{size{G)) and hence, runs in polynomial time for 
a fixed number of vertices. 

Proof. The first step in the algorithm takes time 0{\a{g)2'-'''^^'^°^'^^)poly{svze{G)) since a separation 
oracle for Vm,n can be constructed in polynomial time in the size of G and Kannan's algorithm takes 
size(G)) and we iterate 2"''°s" times. Since h\{g) is polynomially bounded in the size of 
the input, the time complexity of the algorithm is 0(2'^("'°s")poly(size(G)). □ 

3.5 The Algorithm 

We now summarise the results that we obtained in the previous section to obtain an algorithm for 
computing the rank of a divisor. 

Algorithm 5. 1. If deg(D) < then, output r{D) — —1. 

2. If g < deg{D) < 2g - 2 then, set D' = K - D and compute r{D'). Output r{D) = r{D') + 
deg{D)-{g-l). 

3. Ifdeg{D) > 2g - 2, then r{D) = deg{D) - g. 

4- ^/O < deg{D) < g — 1, then we invoke Algorithm^ to compute r{D). 

Remark 5. As described in 11], the notion of rank of a divisor can also be defined for an arbitrary 
sublattice of the root lattice A^. In such a general setting, we do not know if rank can computed in 
polynomial time even when the dimension of the lattice is fixed. In our algorithm we crucially exploit 
our knowledge of the extremal points and the problem with handling the general case is that we do 
not have an explicit description of the extremal points as we have in the case of Laplacian lattices. 
As a consequence, we do not know how to find the extremal points in polynomial time even when the 
dimension is fixed. 

We will end this section by determining the vertices of the polytope Vm^n- 

Lemma 15. The vertices ofVj^-^ ji^ are of the form Wij = Riti — R2tj for i j where to, . . . ,tn are 
the vertices of A . 

Proof. Observe that Wi,j is contained in T^a^.^j for all pairs i, j from to n. We now show that 
is contained in the convex hull of Wi^j where i, j vary from to n. Let p be a point in Vb.i,R2 by 
definition it can written as Ri X]r=o ~ ^2 X]"=o '^i^i with > 0, > for i from to n and 
YIi^o^i = Sr=o^« = 1- We let lij = Xi.aj and write p = Z^-'^q Z)j=o ■ "^^"^ verify that 

J2i j = 1 and £ij > 0. This shows that Vr^^r^ is contained in the convex hull of {wi,j}i^j. 

We now show that Wi^i are not vertices of Vr^^r^ since Wi^i is contained in A(0,i?i — R2) and 
A(0,i?i — R2) is contained in A{0,Ri) if i?i > R2 and is contained in A(0,i?2) otherwise. To 
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conclude the proof of the lemma, it suffices to show that Wij is a vertex if z ^ j. To this end, we 
consider the linear function fij — Xi — Xj and note that Wij is the unique maximum of among all 
points in the set {wij}ij. □ 



4 A duality theorem and its generalisations 

We obtain another proof of the fact that testing if a divisor D is effective or not i.e., testing if r{D) > 
is contained in NP O co — NP. Although, this fact is already implicit in the algorithm of Tardos, our 
proof motivates a generalisation of the duality theorem which will see in the next subsection. 

An interesting characterisation of divisors of negative rank (i.e., rank equal to minus one) is the 
following theorem first established in [5]: 

Theorem 8. A point D has rank minus one if and only if it dominates a point in Ext{LQ). 

Theorem |8] was further generalised to arbitrary full dimensional sublattices of in [1 . We will 
state a "projected" version of the theorem here. Let L be a full-dimensional sublattice of Hq. For 
any real number t > 0, define the arrangements At = UggLA(g,t) and Bt = ^Jc^Crit i^(L)^{c,t,) (see 



Subsection 2.5 for the definition). We have the following duality theorem for the arrangements At and 
Bf 

Theorem 9. A DuaUty Theorem for the arrangement of SimplicesJ Let L he a full dimensional 
lattice in Hq for any real number < t < Cov^{L), then the arrangements At and Bcovi^(L)-t Hq 
i.e., int{At) U Bcov^{L)~t = Hq and int{At) n Bcov^(L)-t = 0- 



We now use Theorem 12 and the fact that in the case of Laplacian lattices, the set of CritA(iG) is 
precisely the projection of Ext{La) onto Hq shows that the computational problem of testing whether 
rank of a divisor is non-negative and is in NP C\co — NP. As we shall see the fact that the elements 
of Giit /^{Lq) have a combinatorial interpretation help to obtain an efficient verification procedure for 
r{D) < 0. 

Theorem 10. The problem of deciding if r{D) > is contained in NP n co — NP. 

Proof. Suppose that r{D) > then by definition there is a point q E Lq such that D > q. We can 
assume that q is contained in A{ttq(D),Covi^{Lg)) and hence, q has a description that is polynomial 
in the description of D and Cov/^{Lg). We know that Covi^{Lg) is the density of the graph (see 
^). The verifier then tests if D > q and whether q G Lq i.e., and hence, if q can be written as an 
integral combination of the basis elements of Lq- This test is the lattice membership problem and can 
be performed in polynomial time, see [Hj for more details. 

This shows that the problem of deciding if r{D) > is in NP. Conversely, to show that the 
problem of testing if r{D) < is also contained in NP, we use the structural theorem of Theorem 
[sj By Theorem [s] we know that if r{D) < then there is a point c' in Grit aILg) of the form c' — 

+ q where q G Lq such that c G 1\{'K{){D),Covi^{Lq) — deg{D)/n + 1). Since, c' is contained in 
1\{'Kq{D),Covi^{Lg) — deg{D)/n + 1) we know that c' has a description that is polynomial in the 
description of D and the description of G. In order to verify that c' is of the form c' = + g for 
some TT G Sn+i and q G Lq, the prover also produces tt and q. The verifier then tests if q is in Lq 
by testing if q can be written as an integral combination of a basis of Lq, this test is known as the 
membership testing problem and can be performed in polynomial time and then test if c' — g is equal 
to the projection of {indegT^(vi) — 1, indegTr{v2) ~ 1, • ■ • , indegjrivn+i) — 1) onto Hq. Hence, deciding if 
r{D) < is also contained in NP. □ 



4.1 Generalised Duality 

Reformulating Corollary[2]in a more geometric language, the rank of the divisor is essentially the smallest 
possible radius of the arrangement of polytopes that are up to scaling the unit ball of the ^i-norm and 
centered at points in Cnt^{Lc) that contains ■kq{D). A "duality theorem" similar to Theoremlolfor the 
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-unit ball instead of the simplex A will, under some mild assumptions, have similar implications as 



Theorem 10 for the problem of testing if a divisor of degree g — I has rank at most c for any c between 
— 1 and g — 1. This observation also raises the problem of extending the notion of duality from the 
regular simplex A to a general polytope. 

Definition 13. (IDuality for polytopesj The graph G is said to have duality with respect to the 
polytope V if there exists a star-shaped body V* such that for R > 0, there exists an R' > such that 
a point p is contained in the interior of the arrangement UceCrit/s{Lc)'^ic, R) if and only if it is not 
contained in [Jq^^^V* {q, R' ) . 

Remark 6. Our imposition that V* has to he star-shaped in the definition of duality is somewhat 
arbitrary and is motivated from the fact that in practice V* happens to be star-shaped. 

While we do not know of tools to obtain duality theorems with respect to polytopes in general we will 
now sketch a general strategy to prove duality theorems for simplices. Here are the main ingredients: 

(11) . From the results in [1], we know that the elements of CritA(iG) have a geometric interpretation 
namely as the local maximum of the simplicial distance function A on the lattice Lq- We use the fact 
that the set of local maxima of the simplicial distance function behave "well" under non-singular linear 
maps i.e., if M is a non-singular linear map then M{Crit/\{LG))) is the set of local maxima of M{Lg) 
under the simplicial distance function dM{A)- 

(12) . Suppose that the linear map M fixes Lq and also CritA(iG) as sets then we know that 
CritA(iG) must be the set of local maxima of the simplicial distance function induced by the simplex 
Af(A) on Lq. This motivates a new notion of automorphism of a graph i.e., the group of linear 
transformations that fix both Lq and CritA(iG)- Furthermore, the graph G has duality with respect 
to the simplex M(A) with (M(A))* being -M(A). 

Let us now discuss each item of the general strategy in more detail. 

Lemma 16. Let AI be a non-singular linear transformation on Hq and let L be a full dimensional 
lattice in Hq, then for every point p G Hq, h^ j^{p) — /im(A).m(l)(-^^(p))- 

Proof. We know that there exists a lattice point q (z L such that h/\ i^(p) — d/\{p,q). Hence, p G 
A(0, hA,L{p)) + q- We have: M{p) € M{A{0, hA,Lip))) + M{q) and we have: dM(A){M{p), M{q)) = 
/ia.l(p) and hence, h]^j(^/x).M{L){M{p)) < hA,L{p)- Since, M is non-singular, is well defined 

and we apply the argument for a point in M{L) and the distance function dM{A) to deduce that 

hM(A)M(L){M{p)) > h^xip)- □ 

Another result we use is the behaviour of arrangements of simplices under non-singular linear trans- 
formations. Let I? be a discrete point set of Hq and let As.'D,t be the arrangement of simplices of type 
S on V, more precisely Av.s.t ~ ^qevS{q,t). 

Lemma 17. For any non-singular linear transformation M on Hq, M{Av,s,t) — AM{v),M(S).t ■ Fur- 
thermore, the map M preserves the interior of the arrangements i.e., M{int{Av,s,t)) = int{A(M(v),M{s),t))- 

Proof By definitiom M{Av,s,t) = M(U,ei,5(0, i) + g) = (Uqei,M(^(0, t) + g) = (U,ei,M(5(0, t))) + 
M{q) — AM{T>),M(S),t- Now, suppose p is a point in the interior of the arrangement ^x>,s,t then we know 
that there exists e > such that B{p,e) in the arrangement AT>,s,t- Applying the linear transforma- 
tion M, we know that M{B(jp,t)) is contained in the arrangement M{Axi,s.t) = A(^M{v),M{S),t)- Now, 
we know that M{B{p,e)) contains the ball B{M{p), ^A(M)e) where A(M) = inf{||M(a;)|||| \\\x\\'^ = 
1, X £ Hq} since M is non-singular we know that X{M) > and hence, M{p) is contained in the 
interior of A(M{V).M{S),t) and similarly, if we apply the same argument for a point p' in the inte- 
rior of A(^M{V).M(S),t) then M~^{p) is a point in the interior of Ax>,s,t- Hence, M{int(Ax>,s,t)) = 

/nt(^(jV/(X)),M(5),t))- Q 

Corollary 4. Let M be a non-singular linear transformation on Hq and let L be a full dimensional 
lattice in Hq, then Gova{L) = CovM(h.){M{L)). 

Remark 7. The second part of Lemma \l'}\ is essentially a topological fact and in topological terms we 
use the fact a non-singular linear transformation from Hq to itself is a homeomorphism from Hq to 
itself where Hq is equipped with the Euclidean topology. 
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Remark 8. Note that we did not use our assumption that the set V is a discrete point set anywhere in 
the proof of Lemma 17 and can be generalised to non-discrete sets but since we mainly deal with discrete 
sets we restrict to them. 



Theorem 11. Let M be a non-singular linear transformation on Hq and let L be a full dimensional 
lattice in Hq, the set of local maxima of the function /im(A),M(L) *s precisely the set M{Crit^(L)). 

Proof. Suppose c is an element of Crit/\(L), we know that there exists an e > such that h^x{c) > 
ft.A.L(c') for all c' e B{e,c). Now by Lemma 16 we know that /ij\/(A),j\/(L) (-^^(c)) > ft-M(A),Af(L) (c") 



for all c" G M(B(e,c)). The set M(B(e,c)) is an ellipsoid that contains the baU B{M{c),e^/X{M)) 
where A(M) is inf{||M(a;)||2| \\\x\\2 = 1, x G Hq} and since M is non-singular we have A(M) > 0. 
Hence, M{Crit/\{L)) C Crit m( A) {M{L)). Since, M is a non-singular linear map, we consider the map 
and apply the above argument for a point in CTitM(A){^'I(L)) to deduce that Crit^/(-/x) (Af (L)) C 
M(CritA(i)). □ 

We will now generalise the duality theorem with respect to any full dimensional simplex on Hq with 
centroid at the origin. 

Theorem 12. Let S be a full dimensional simplex on Hq with centroid at the origin and let AL,s,t be 
the arrangement U^gi — S{q,t) and let BL,s,f be the arrangement Ucec'rits(L)S{c,t). Fix a real number 
<t < Cova{L), the arrangements A L,s,t and Bl,s, Govs {L)-t tile Hq i.e., AL,s,t^BL^s,Covs{L)-t = Hq 
and int{AL,s,t) n BL,s,Covs{L)-t = 0- 

Proof. Since 5 is a full-dimensional simplex with centroid at the origin we know that there is a non- 
singular map M that takes S to the simplex A. Furthermore, by Corollary [4| we know that Cov-s{L) = 
Cova{M{L)) and hence < t < Cov^{M{L)) and we can apply Theoremlolto the lattice M{L), a full 
dimensional lattice in Hq and with the simplex A, that a point p G Hq is contained exclusively in either 
exclusively the interior of the arrangement AM(L),A,t or the arrangement BM(L),A,Cov^(M(L))-t- Now 



note that M{L) and CritA(Af (L)) are discrete sets and apply Lemma 17 with the linear transformation 
M"^ Hence, we know that mi(^i,,s,t)Ui3i^S^C'oi;s(i)-t = Hq and int(AL,s,t)f^BL^s,Covs{L)^t = since, 
if there is a point p € Hq in int{AL,s,t) H BL^s,Covs(L)-t then M~^{p) is contained in int{AM(L),A,t) H 
B(M{L),A,Cov^{M{L))-t) which is a contradiction. □ 



4.2 The Automorphism Group of Critp(L) 

We now discuss the second item of the general strategy, the notion of automorphisms: 

Definition 14. (^Automorphism of a lattice^ For a lattice L, a non-singular linear map M on the 
linear space spanned by the elements of L is called an automorphism of L if it induces a bijection on L. 

Here is a concrete characterisation of the automorphism group of a lattice: 

Lemma 18. Let B be an arbitrary basis of a n-dimensional lattice L written as a matrix with the basis 
elements row-wise. A linear transformation M is an automorphism of the lattice L if and only if it 
M{B) = T-B where T e SL{n,Z). 

The set of automorphisms of a lattice has a natural group structure under composition (if you think 
of the linear maps as matrices then the group operation is matrix multiplication) and this group is 
called the automorphism group of L. 

Definition 15. (Automorphism of Crit a (-^) J Let L be a sublattice of Hq, a non-singular linear 
map M on the linear space spanned by the elements of L is called an automorphism of Crit a (L) if it 
induces a bijection on CritA^L). 

The set of automorphisms of Crit a (L) has a natural group structure under composition (if you think 
of the linear maps as matrices then the group operation is matrix multiplication). The group gives rise 
to a graph invariant that we call the critical automorphism group of the graph. 
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Remark 9. Note that the notion of a critical automorphism group can he developed more generally, 
for any distance function induced by a convex body V on a lattice L, we can define the group of trans- 
formations that fix both L and Critp{L) as sets. 

A convenient way to think of the automorphisms of CritA(iG) is to consider Grit /:^{Lq) / Lq as a 
subset of the torus T^j = Hq/Lq and observe the automorphisms of Lq induce an bijection of the 
torus Tc to itself; among these bijections the automorphisms of CritA(iG) are those that take the set 
Ciit/\{La)/ Lq to itself. Let us now give some concrete examples of CritA(iG): 



1. For a regular graph G, we have CritA(iG) — ~CritA(iG)- Hence, the linear transformation 
M{x) = —X is an example of an automorphism of CritA(iG)- 

2. The permutation map induced by every automorphism of also is indeed an automorphism of 
CritA(iG); but since permutation maps fix the simplex A, these are not useful to obtain simplices 
other than A for which duality holds. Note that these automorphism take one class in Ciit ^{Lq) 
to another. 

3. Every c G CritA(iG) where c ~ c-^ + q and q G Lq. We can write each c^r as a rational 
linear combination of {6i,...,6„}. Define the height of c^r as the least common multiple of 
the denominators of these rational coefficients in their reduced form and define the height of 
CritA(iG) as least common multiple of the ratio of the denominators of the height of each Ct^. A 
linear transformation M{hi) — bi+H{J2j^i o^jbj) is an automorphism of Xg where H is a. multiple 
of the height of CritA(iG) and a^s are arbitrary integers is an automorphism of Crit /x{Lg) / Lq . 
These automorphisms takes an equivalence class Crit /\{Lq) / Lq to itself. 

The usefulness of the critical automorphism group arises from the following partial characterization 
of simplices (with centroid at the origin) for which duality holds: 

Theorem 13. Let G be an undirected connected graph, duality holds with respect to a simplex S with 
centroid at the origin if there is a linear map M taking A to S that belongs to the automorphism group 
ofCrit^{LQ). 

Remark 10. We do not know if the converse of Theorem \13\ also holds, the main difficulty in proving 
the converse is to show that if S has duality then Crit/\(LQ) is the precisely the set of local maxima 
under the distance function ds ■ 



The fact that an automorphism M of CritA(iG) fixes CritA(-/jG) combined with the Theorem 12 
means that the graph G has a duality with respect to the simplex M(A) with the dual (Af(A))* being 
-M(A). 

4.3 The case of complete graphs 

We will study the critical automorphism group of the complete graph. Let 6o, . . . ,6„ be the rows of 
the Laplacian matrix of the complete graph i^n+i. The extremal points are permutations of the point 
(—1,0, 1, . . . ,7i) of the coordinates translated by a lattice point i.e., a point in the Laplacian lattice 
of Kn+i. Since, the points in CritA(iG) are orthogonal projections of the extremal points onto the 
hyperplane Hq = (1, . . . , l)-*-, they are points of the form c^^ + q where = 7r((n/2, n/2 — 1, n/2 — 
2, . . . , —n/2)), where tt is a permutation matrix. Now, since we now that up to equivalence modulo Lq 
there are only n! class of CritA(iG) and that permutations such that 7r(n) = n form representatives of 
the set CvitAiLQ) / Lq. We now consider permutations such that 7r(n) = n and for such a permutation, 
we can rewrite in the basis of {5o, . . . , 6„} as 

Ctt = • (8) 

n 

Hence, the height of CritA(-ftrn+i) is n. We can use the exploit Equation ([s]) to obtain a better 
formula for rank. We will now construct critical automorphisms of Kn+i: 
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Lemma 19. Every map of the form hi = 67r(i) where tt is a permutation is an element of the critical 
automorphism group. 

More generally, we have the following result: 

Theorem 14. A map M that takes hi to h^^^^i^ + H ■ qi for i from to n ~ 1 where n is a permu- 
tation and H is an integer that divides n, the height of Kn+i and qi is contained in the sublattice 
spanned by 67^(0), • • • , ^77(4-1) i-^-j Qi = Sj=o '^ij^-^U) where aijS are integers is an element of the critical 
automorphism group of K^+i- 

Proof. First, we show that the map takes Lx^j_i to itself as follows; we write for * from to n as 
an integral combination of M(67r(o))i ■ • ■ : -^'^(^7r(n)) ^ follows: 6,^(0) = M{bi) and for i from 1 to n — 1 
we have: 



= M{bi) ^okMibj) (9) 

j=0 fc=0 

Hence, M takes a basis of Lk,^^^ to another basis of Lk„j_i and hence takes Lk„+i to itself. 

Now, to show that M takes CritA(ii<-„+i) to itself we consider M[ct^) and form Equation (jsj) we 
have: 

M{c„) = — — — — ^ ^j. Now smce H divides n we know that 

M{ca) is mapped to c^on + q' for some lattice point q' and since M takes Kn+i to itself it takes 
M^c^+Lk^^i) to CcroTT+iifn+i- Furthermore, since tt and a are permutations, M induces a permutation 
on equivalence classes in CritA(iif„+i)/i_fs'„+i represented by with 7r(n) = n. These arguments show 
M that takes CritA(-£'_R:„+i) to itself and hence M is an element of the critical automorphism group of 

Kn+l- " □ 

Remark 11. Note that Theorem \14\ provides a construction of elements of the critical automorphism 
group of Kn+l, but the construction is not exhaustive in the sense, there are other elements of the 



critical automorphism group of Kn+i that are not constructed by Theorem I4 



4.4 Concluding remarks 

A natural outgrowth of our work would be to obtain a singly exponential algorithm for computing the 
rank and furthermore, to resolve the complexity of computing the rank i.e., to determine whether it is 
NP-hard or it is computatble in polynomial time. The bottleneck in our algorithm is the enumeration 
of all permutations. We believe that another interesting direction would be obtain duality theorems for 
polytopes generalising the results that we obtained on simplices. 

Acknowledgements: The author thanks Omid Amini, Khaled Elbasionni, Jan van den Heuvel 
and Amr Elmasary for the stimulating discussions they had with him on the topic of the paper. 
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